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Abstract. This study evaluates the spatial distributions and seasonal cycles in upper 
tropospheric ozone (pressure range 200-500 hPa) from low to high latitudes (60S to 60N) 
derived from the satellite retrieval method called “Cloud Slicing.” Cloud Slicing is a 
unique technique for determining ozone profile information in the troposphere by 
combining co-located measurements of cloud-top pressure and above-cloud column 
ozone. For upper tropospheric ozone, co-located measurements of Nimbus 7 Total Ozone 
Mapping Spectrometer (TOMS) above-cloud column ozone and Nimbus 7 Temperature 
Humidity Infrared Radiometer (THIR) cloud-top pressure during 1979-1984 were 
incorporated. In the tropics, upper tropospheric ozone shows year-round enhancement in 
the Atlantic region and evidence of a possible semiannual variability. Upper tropospheric 
ozone outside the tropics shows greatest abundance in winter and spring seasons in both 
hemispheres with largest seasonal and largest amounts in the NH. These characteristics 
are similar to lower stratospheric ozone. Comparisons of upper tropospheric column 
ozone with both stratospheric ozone and a proxy of lower stratospheric air mass (i.,e., 
tropopause pressure) from National Centers for Environmental Prediction (NCEP) 
suggest that stratosphere-troposphere exchange (STE) may be a significant source for the 
seasonal variability of upper tropospheric ozone almost everywhere between 60S and 
60N except in low latitudes around 10S to 25N where other sources (e.g., tropospheric 
transport, biomass burning, aerosol effects, lightning, etc.) may have a greater role. 

1. Introduction. 

The first assessment of tropospheric column ozone (TCO) derived from satellite 
measurements was based on the residual method of Fishman et al. [1990] which 
subtracted SAGE stratospheric column ozone (SCO) from TOMS total col umn ozone to 
determine TCO. Although it was demonstrated that this method could retrieve 
information of TCO distributions extending from the tropics to middle latitudes, 
essentially only a climatology of TCO was obtainable because of the sparse solar 
occultation measurements of stratospheric ozone profiles available from SAGE (around 1 
to 2 local profile measurements per month). In addition it was noted by Fishman et al. 
[1990] that this method had further difficulties caused by subtracting two large column 
ozone measurements from two independent instruments (i.e., TOMS and SAGE) which 
were not inter-calibrated. 

Subsequently, a number of algorithms based on tropospheric ozone residual methods 
were developed with significant improvements in characterizing the seasonal, interannual, 
and spatial variations of TCO in the tropics [e.g., Jiang and Yung, 1996; Kim and 
Newchurch, 1996; Hudson and Thompson, 1998; Ziemke et al., 1998; Thompson and 



Hudson, 1999; Fishman and Balok, 1999; Kim et al, 2001; Ziemke et al., 2001]. The 
TCO derived from these techniques have yielded many important scientific results 
including characterization of ozone enhancement in the Atlantic during southern spring 
[first identified by Fishman et al., 1987], evidence of an El Nino signal [Chandra et al., 
1998, Thompson et al., 2001], solar cycle forcing [Chandra et al., 1999], characterization 
of seasonal and interannual variabilities [Ziemke and Chandra, 1999], identification of 
lightning-generated TCO distributions [Martin et al., 2001], and a tropical Atlantic 
Paradox [Thompson et al., 2001]. Additional studies have shown that the September- 
October Atlantic maximum and/or the El Nino increase of tropospheric ozone over 
Indonesia in 1997 were well captured in 3D photochemical transport models [e.g., Moxim 
and Levy, 2000; Sudo et al., 2001; Bey et al., 2001; Martin et al., 2002; Chandra et al., 
2002a]. 

While much has been learned about tropical TCO from various satellite retrieval 
methods, the extension of measurements and the characterization of TCO outside the 
tropics has had limited success. The residual method of Fishman et al. [1990] was more 
recently applied by replacing SAGE measurements with solar backscatter ultraviolet 
(SBUV) to determine SCO [Vukovich et al., 1996; Fishman and Balok, 1999], Unlike 
SAGE the SBUV profiles provide daily measurements of SCO over much of the globe. 
However SCO derived from SBUV has limited success because of poor vertical 
resolution of ozone profile information in the vicinity of the tropopause. Fishman and 
Balok [1999] provided an improvement to this TOMS/SBUV residual technique by 
applying a normalization to the SBUV profile measurements using ozonesonde 
climatology [Logan, 1999]. 

The Cloud Slicing technique [Ziemke et al., 2001] provides a method for deriving ozone 
in the upper troposphere over much of the Earth wherever sufficient clouds exist and co- 
located measurements of above-cloud column ozone and cloud-top pressures are 
available. The Nimbus 7 TOMS instrument determines total column ozone in the 
atmosphere from measuring backscattered UV radiances at six wavelengths varying from 
3 13 nm to 380 nm. , Optically thick water vapor clouds in the troposphere are opaque at 
these UV wavelengths, and as a result the column ozone amount actually measured by 
TOMS in the presence of such clouds is “above-cloud” column ozone. Cloud Slicing 
takes advantage of the opaque nature of these thick clouds and above-cloud column ozone 
measurements to estimate ozone abundance in the troposphere. With this method, above- 
cloud column ozone is plotted versus co-located cloud-top pressure for a chosen 
geographical region and for a pre-selected pressure band. The mean slope of the 
distribution directly yields mean ozone volume mixing ratio. 

The current study derives upper tropospheric ozone volume mixing ratio in the pressure 
range 200-500 hPa between latitudes 60S to 60N and characterizes both spatial patterns 
and seasonal variabilities present in the data. Although not the main objective, we also 
investigate possible implications for stratosphere-troposphere exchange (STE) by 
comparing spatial distributions and seasonal cycles in upper tropospheric ozone with 
lower stratospheric ozone from UARS Halogen Occultation Experiment (HALOE), 



TOMS total column ozone measurements, and lower stratospheric air mass derived from 
NCEP re-analyses. Ozone in the upper troposphere and lower stratosphere exhibits 
lifetimes of months to years, respectively. In these regions ozone is a chemical tracer of 
air mass motion and distribution. One may expect that upper tropospheric ozone, lower 
stratospheric ozone, and lower stratospheric air mass will show similar seasonal cycles 
with the occurrence of STE. 

The following sections begin with a description of the data and an overview of the Cloud 
Slicing method, followed by ozonesonde validation of measurements, characterization of 
seasonal cycles and spatial distributions, correlations of seasonal cycles, and finally a 
Summary. 

2. Nimbus 7 Data and Cloud Slicing. 

Nimbus 7 Temperature Humidity and Infrared Radiometer (THIR) cloud-top pressure and 
Nimbus 7 TOMS version 7 column ozone level-2 footprint measurements were combined 
to determine upper tropospheric ozone from the method of Cloud Slicing [Ziemke et al., 
2001], The THIR instrument provides thermal emission measurements in a 2 pm- wide 
channel centered about 1 1.5 pm to determine cloud-top temperature. THIR cloud-top 
pressure was derived from temperature using National Centers for Environmental 
Prediction (NCEP) analyses. Monthly ensemble averages from Cloud Slicing were 
derived from daily values from January 1979 through October 1984 between latitudes 
60S and 60N (5 degree by 5 degree resolution). The pressure band for Cloud Slicing was 
chosen following several trial data runs to be 200-500 hPa to optimize both number of 
measurements and signal-to-noise. In tropical latitudes the tropopause lies around 100 
hPa year-round while in middle and high latitudes the tropopause may lie at 300 hPa or 
greater pressure. Cloud Slicing measurements represent mean ozone volume mixin g ratio 
in the upper troposphere lying between 500 hPa and highest cloud tops limited to 200 hPa 
in THIR cloud pressure data. In annual cycles the derived distributions of ozone may to 
first order be interpreted as estimates of mean ozone volume mixing ratio lying between 
500 hPa and the tropopause under the premise that ozone is relatively well mixed (i.e., 
constant mixing ratio) on long time scales in the upper troposphere up to tropopause 
level. Further details for the Cloud Slicing algorithm are discussed in Appendix A. 

Figure 1 shows scatter plots of TOMS above-cloud column ozone versus co-located 
THIR cloud-top pressure for a grid point located near Wallops, Virginia (38N, 75 W). 

The top three frames are three days in July 1984 and the bottom three frames are days in 
August 1 984. The numbers shown in each scatter plot are the derived ozone volume 
mixing ratio and 2-sigma uncertainties (in parentheses) in units ppmv. The scatter in 
these distributions is caused by errors in TOMS above-cloud column measurements and 
also errors caused by some amount of stratospheric ozone variability over the selected 5 
degree by 5 degree region (see Appendix). Some additional scatter in the plots in Figure 
1 is caused by errors in converting THIR cloud-top temperature measurements to pressure 
using NCEP analyses. Future satellite instruments with greater spectral resolution may 
enable the determination of cloud-top pressure directly instead of having to convert 



measured cloud-top temperature to pressure from meteorological analyses. Methods 
capable of deriving tropopause directly from backscattered UV and visible wavelength 
radiances include the molecular oxygen dimer (02-02) method [S. ****, personal 
communication, 2002] and the rotational Raman scattering “Ring Effect” [Joiner and 
Bhartia, 1995]. 

Figure 2 shows the average number of Cloud Slicing measurements per month in this 
investigation. Several regions have average monthly sampling rates of 3 to 4 values per 
month which is small, yet still comparable to most ground-based ozonesonde stations 
(section 3). There are two regions in Figure 2 where there are no measurements at all 
from Nimbus 7 Cloud Slicing because of an insufficient number of high clouds or high 
reflectivity scenes: (1) the south Pacific oceanic region west of South America, and (2) 
the south Atlantic oceanic region lying between South America and Africa. 

3. Validation with Ozonesonde Measurements. 

Time series of upper tropospheric ozone were compared with similar time series at 
several ground-based World Ozone and Ultraviolet Radiation Data Center (WOUDC) 
stations. Figure 3 shows time series comparisons of upper tropospheric ozone at several 
ground station locations. It is noted that both data sources were constructed as monthly 
ensembles which do not represent true monthly means. This is because measurements for 
both sources are scarce with some monthly ensembles having as little as only one day of 
measurements. Despite data sampling problems, upper tropospheric ozone time series in 
Figure 3 indicate general agreement in mean amounts and seasonal variability with 
greatest abundance around spring months and minimum amounts in Aut umn in both 
hemispheres. Figure 4 shows seasonal cycle comparisons of the data plotted in Figure 3, 
including comparisons at several additional ground-station locations. The station 
locations in Figure 4 begin in the upper left-most frame for Goose, Canada (53N, 60W) 
and progress with decreasing latitude from left-to-right and downward by row. Most all 
of the stations are in the NH extratropics (i.e., Europe, Japan, U.S.) and show largest 
seasonal variations and mean column amounts in mid-to-high latitudes. The SH stations 
(last two frames) show a reversal of the seasonal cycle variability with largest ozone 
around SH spring months. 

4. Characterization of Seasonal Cycles and Spatial Distributions. 

The primary purpose of this investigation is to extend Cloud Slicing to middle and high 
latitudes and characterize observed seasonal cycles and spatial distributions of upper 
tropospheric ozone. Figure 5 shows spatial distributions of 6-year averaged 3-month 
seasonal means (indicated) of upper tropospheric ozone volume mixing ratio. Largest 
ozone occurs in winter and spring months in both hemispheres with greater abundance in 
the NH. The most significant amounts occur in mid-to-high latitudes. In the NH 
midlatitudes in winter and spring months there is evidence of partial poleward entrapment 
of upper tropospheric ozone caused by the upper tropospheric wind jets which can act as 
a partial barrier for meridional transport of air mass in the upper troposphere. In the NH 



in winter and spring there are large meridional gradients and enhancements in ozone 
along wind-driven storm track regions such as the eastern U.S. and the Asian continent 
extending into the eastern Pacific. In the SH in winter and spring seasons the horizontal 
gradients and mean amounts in ozone are weaker than winter and spring seasons in the 
NH. There is also a persistent feature present in each season in the SH subtropics which 
indicates an enhancement of ozone near South America, Africa, and Australia. In the 
tropics, enhancement of upper tropospheric ozone is present each season in the Atlantic 
region from South America to Africa compared to the eastern and western Pacific 
regions. 

The seasonal and spatial characteristics of upper tropospheric ozone in Figure 5 resemble 
that of total column ozone (Figure 6). We note that the near-global similarity of upper 
tropospheric ozone and total ozone is not an artifact of the Cloud Slicing algorithm which 
assumes that stratospheric column ozone is spatially invariant over selected 5 degree by 5 
degree regions (i.e., all above-cloud column ozone variability over the region is 
tropospheric in origin). Contamination from stratospheric ozone variability within the 
algorithm is limited to intense dynamical regions, particularly the tropospheric wind jet 
regions which are limited to narrow meandering bands in the subtropics and middle 
latitudes. The Cloud Slicing algorithm removes outlier data associated with the wind jet 
regions by applying statistical and column ozone filtering (see Appendix A). 

The most coherent similarities between upper tropospheric ozone in Figure 5 and total 
ozone in Figure 6 occur in the NH in winter and spring months. Most contribution to 
total column ozone lies in the mid-to-lower stratosphere. Because ozone is generally a 
long-lived tracer of air mass in both the upper troposphere and lower stratosphere 
(lifetimes of months to years, respectively), the patterns in Figures 5 and 6 suggest a 
possible influence from STE (discussed in section 5). Figures 7 and 8 show zonal mean 
seasonal cycles from the data plotted in Figures 5 and 6. It is noted for the upper 
tropospheric ozone data in Figure 7 that these are not true zonal means since in several 
months there are regions without enough suitable clouds for applying the Cloud Slicing 
method (e.g., see Figure 5 and discussion). Nevertheless, in Figures 7 and 8 both upper 
tropospheric ozone and total column ozone show similar seasonal variation and latitude 
dependence. In the tropics seasonal cycles are different in these quantities, and both 
temporal and meridional variations appear weak. Outside the tropics, mean amounts and 
meridional gradients are greatest around spring months in both hemispheres. 

Figure 9 shows zonal mean upper tropospheric ozone for January 1979-October 1984. As 
in Figure 7, variabilities in upper tropospheric ozone in Figure 9 show largest latitudinal 
gradients in the extratropics of both hemispheres with peak amounts around spring 
months. It is noted that the springtime maximum in upper tropospheric ozone volume 
mixing ratio in NH middle and high latitudes agrees with the NH ozonesonde assessment 
of Logan [1999] (e.g., her Figure 27 for 300 hPa in mid-to-high latitudes and 200 hPa in 
the subtropics). Figure 9 indicates a possible semiannual component in the tropics. Peak- 
to-peak amplitudes of this signal are weak, around 2-4 ppmv. In some years (e.g., 1980- 
1983) there is evidence of a semiannual oscillation while in other years there is not. A 



semiannual component in the tropics may indicate a dynamical and/or photochemical 
forcing related to the semiannual zenith crossing of the sun in the tropics. Comparisons 
of the semiannual and annual variations in total column ozone and upper tropospheric 
ozone are shown in Figure 10. The annual and semiannual amplitudes were derived by 
applying basic Fourier analysis. Both annual and semiannual components in total ozone 
and upper tropospheric ozone in Figure 10 compare well in a global sense with annual 
cycle variations exhibiting larger amplitudes and larger latitudinal gradients. 

5. Correlations of Seasonal Cycles. 

It has been shown that total ozone and upper tropospheric ozone indicate similar spatial 
and seasonal cycle variabilities extending from low to high latitudes. A plausible 
explanation for these coherent temporal and spatial relationships may be STE processes 
in the atmosphere. It is well known that ozone in the upper troposphere and lower 
stratosphere is a relatively long-lived chemical tracer (e-folding lifetimes of months to 
years, respectively) so that even monthly measurements in this study provide valuable 
insight regarding possible occurrence of STE. Recent studies suggest that tropopause 
folding and turbulent convection are significant sources for transporting stratospheric air 
mass into the troposphere [e.g., Langford et al., 1996; Beekmann et al., 1997; Cho et al., 
1999; Eisele et al, 1999; Meloen et al., 2001; Bertin et al., 2001; Fujiwara and 
Takahashi, 2001]. The global effects of STE on tropospheric ozone and its seasonal 
variability may be significant. The 3D model study by Lelieveld and Dentener [2000] 
indicated that upwards of 80% of tropospheric ozone may be ascribed to STE in 
subtropical-to-high latitudes. 

Appenzeller et al. [1996, and references therein] showed that air mass lying in the 
lowermost stratosphere is fundamental to global STE processes. This study investigates 
potential STE on observed upper tropospheric ozone seasonal cycles by comparing with 
both lower stratospheric ozone and lower stratospheric air mass. Air mass abundance in 
the lowermost stratosphere is parameterized from tropopause pressure (see Appendix B). 
For lower stratospheric ozone we include measurements from the UARS HALOE 
instrument. HALOE is a solar occuitation instrument which provides profile inf o rmat ion 
of stratospheric ozone from the top of the atmosphere down to the tropopause [. Bruhl et 
al., 1996]. Mean ozone volume mixing ratio from HALOE was calculated for both the 
46-100 hPa and 46-215 hPa pressure bands (i.e., together encompassing the lower 
stratosphere extending from tropics to high latitudes) for the 1993-2000 post-Pinatubo 
time period. Although for a time period later than 1979-1984 and with sparse temporal 
and spatial coverage as an occuitation instrument, the combined 8 years of HALOE 
measurements provide retrievals of lower stratospheric column ozone adequate to 
develop seasonal cycles for this comparison. Figure 1 1 shows seasonal cycles in HALOE 
zonal mean ozone volume mixing ratio for the 46-100 hPa pressure band evaluated for 
1993-2000. For latitudes less than around 25-30 degrees the 46-100 hPa pressure band 
represents the lower stratosphere. At higher latitudes the pressure band is more 
representative of the mid-to-lower stratosphere. Seasonal variabilities of HALOE ozone 
in Figure 1 1 are comparable to seasonal variabilities in upper tropospheric ozone in mid- 



to-high latitudes (compare Figure 1 1 with Figure 7). 

Figure 12 shows temporal correlations as a function of latitude between upper 
tropospheric ozone and HALOE (stars, long-dashed curves), TOMS total column ozone 
(bold), and NCEP tropopause pressure (dotted). Correlation calculations with 46-215 hPa 
HALOE measurements were not plotted in Figure 12 equatorward of 35 degrees latitude 
because 215 hPa exceeds tropopause pressure in these lower latitudes in winter and 
spring months. The correlations in Figure 12 are based on 12 months of zonal mean data 
at each latitude and therefore prescribe only basic agreement or disagreement between 
seasonal cycles. Correlation amplitudes exceeding 0.58 pass 95% confidence level (i.e., 
5% significance level). Seasonal cycles for upper tropospheric ozone, total ozone, and 
NCEP tropopause pressure were determined from 1979-1984 data. Seasonal cycles for 
HALOE lower stratospheric column ozone were evaluated for the 1993-2000 (i.e., post- 
Pinatubo) time period. 

If STE is the dominant source of ozone seasonal variability in the upper troposphere, then 
all four curves in Figure 12 would conceivably have correlations near +1 at latitudes 
where this persists. In subtropical-to-high latitudes, NCEP tropopause pressure and both 
HALOE lower stratospheric ozone and TOMS total ozone measurements all indeed show 
correlations near +1 . However this is not the case in the low latitudes in Figure 12 
extending from around 10S to 25N. In these lower latitudes the correlations with NCEP 
tropopause pressure become negative in the NH tropics and correlations with HALOE 
and TOMS are greatly reduced to low-to-negative values. The poorest correlations occur 
in the NH tropics and are caused by a seasonal maximum in upper tropospheric ozone 
during NH spring season compared to seasonal maxima around NH Autumn months in 
NCEP tropopause pressure, TOMS total ozone, and HALOE lower stratospheric ozone 
(note seasonal cycle differences for mean ozone volume mixing ratio in Figures 7 and 1 1 
in the NH tropics). Several potential sources of upper tropospheric ozone could explain 
the seasonal cycles for latitudes 10S to 25N including biomass burning, biogenic 
emissions, aerosol influenced photochemistry, tropospheric transport, and cloud lightning 
[e.g., Moxim and Levy, 2000; Lelieveld and Dentener, 2000; Bey et al., 2001; Thompson 
et al., 2001; Martin et al., 2001, 2002]. 

6. Summary. 

This investigation has examined spatial distributions and seasonal cycles in upper 
tropospheric ozone (pressure range 200-500 hPa) extending from low to high latitudes 
(60S to 60N) by applying the satellite retrieval method known as “Cloud Slicing“. This is 
a new technique for determining ozone profile information in the troposphere by 
combining co-located measurements of cloud-top pressure and above-cloud column 
ozone. Previous investigations have had limited ability in measuring tropospheric ozone 
beyond tropical latitudes from satellite data. Comparisons between upper tropospheric 
ozone volume mixing ratio derived from Cloud Slicing and ground-based ozonesondes 
for 60S to 60N compare well in monthly ensemble time series and seasonal cycles despite 
sparse measurements available from either data source. These comparisons indicate the 



efficacy of this method even as applied to the relatively coarse field-of-view scanning 
measurements from Nimbus 7 TOMS (50 km by 50 km at Nadir, 100 km by 100 km on 
average). Future satellite instruments with smaller footprint size and greater spectral 
resolution than TOMS may provide significant improvements to the Cloud Slicing 
technique, with larger number of measurements and application to constituents other than 
ozone. 

Upper tropospheric ozone indicates a possible semiannual component variability in the 
tropics. A semiannual oscillation appears in some years but not in other years which may 
be a result of poor signal-to-noise with the Nimbus 7 measurements. The existence of a 
semiannual component in the tropics may indicate a dynamical and/or photochemical 
forcing related to the twice-per-year zenith crossing of the sun in tropical latitudes. 

Spatial patterns in upper tropospheric ozone are similar to total column ozone, 
particularly in the NH during winter-spring. Given long lifetimes (months to years, 
respectively) for ozone in the upper troposphere and lower stratosphere, ozone in these 
regions of the atmosphere is a tracer of air mass movement and distribution. The spatial 
patterns and seasonal cycles in mid-to-high latitudes in both hemispheres are similar for 
TOMS total ozone, lower stratospheric ozone from UARS HALOE, NCEP lower 
stratospheric air mass, and upper tropospheric ozone from TOMS Cloud Slicing. The 
mean amounts and seasonal variabilities in upper tropospheric ozone in the extratropics 
are both larger in the NH compared to the SH. These observations along with previous 
2D and 3D model results indicate a presence of STE in the extratropics. Our analyses 
suggest support for the concept of a coupling of the effects of air mass in the lowermost 
stratosphere (i.e., “middle world“) and the meridional circulation which in model results 
produces largest STE in the extratropics of both hemispheres around spring, minimnm 
STE in autumn, and with greater occurrence in the NH 


Comparisons between lower stratospheric air mass, stratospheric ozone, and upper 
tropospheric ozone imply that STE does not appear to be a major source for the seasonal 
variability of ozone in the upper troposphere in low tropical latitudes extending over a 
broad region from around 10S to 25N. In these latitudes there is poor temporal 
correlations in seasonal cycles between troposphere and stratosphere, especially in NH 
tropical latitudes were largest amounts of upper tropospheric ozone from Cloud Slicing 
occur in spring months. There are several sources for upper tropospheric ozone that 
could explain the seasonal cycles in low latitudes. These sources include biomass 
burning, biogenic emissions, aerosol influenced photochemistry, tropospheric transport, 
cloud lightning, etc. (beyond the extent of this study). 

Appendix A: Cloud Slicing Algorithm. 

Cloud Slicing was applied to daily ensemble measurements over 5 degree x 5 degree 
binned regions, and then averaged within each month for the January 1979-October 1984 
time period. The chosen latitude coverage of co-located TOMS and THIR footprint 
measurements extends from 60S to 60N. This reduces difficulties of persistent snow and 



ice in the polar regions and also eliminates problems with polar night conditions where 
there are no TOMS backscattered UV ozone measurements. For the Cloud Slicing 
technique, footprint measurements (50 km x 50 km at Nadir, around 100 km by 100 km 
on average) of above-cloud column ozone from Nimbus 7 TOMS are plotted versus co- 
located measurements of Nimbus 7 THIR cloud-top pressure over a pre-selected pressure 
band (200-500 hPa in this study). The mean slope of the distribution then directly yields 
mean volume mixing ratio for ozone within that pressure band. Specifically, mean 
volume mixing ratio (X, in ppmv) is determined from 

X = 1 .27 AQ/AP, (Al) 

where Q is above-cloud column ozone (in DU), and P is cloud-top pressure (in hPa). Eck 
et al. [1987] showed for Nimbus 7 TOMS that footprint scenes with reflectivity greater 
than 0.4 are almost always 100% cloud-filled (i.e., 100% cloud fraction scenes). Only 
scenes with reflectivity greater than 0.44 were selected in this study to eliminate most 
partially cloudy footprint scenes. This reflectivity filtering is a relaxation from the more 
stringent 0.6 value of reflectivity used by Ziemke et al. [2001] and provides a greater 
proportion of scenes for applying Cloud Slicing. Another important parameter is the 
minimum number of daily co-located measurements of ozone and cloud pressure within 
each 5 degree x 5 degree region. This number was selected (for adequate statistics) in the 
present study to be at least 20, and as in Ziemke et al. [2001], only computed miying ratio 
values greater than the statistical 2-sigma value were retained. 

A fundamental assumption with Cloud Slicing is that stratospheric column ozone is 
invariant over the selected region (i.e., all above-cloud column ozone changes over the 
region are tropospheric in origin). This condition breaks down in and near strong 

dynamical regionsvparticularly the troposphericwind jet regiohs- where stratospheric - 

column ozone may exhibit large zonal and meridional variability. An additional filter is 
applied to identify these regions: Following analysis for each daily 5 degree x 5 degree 
gridded region, only horizontal gradients in above-cloud column ozone (|AQ|) less than 
50 DU were retained in the data base. This number was chosen subjectively following 
several trial runs with Nimbus 7 Cloud Slicing. It is likely that there is some 
stratospheric contamination in the retrievals, but these effects are expected to be limited 
mostly to narrow regions associated with the tropospheric wind jets. The consistency 
with ozonesonde measurements (section 3) and lack of obvious outliers in these regions 
suggest that the algorithm does not produce a significant amount of stratospheric ozone 
gradient contamination in the retrievals. It is noted that it is also possible to identify the 
wind jet regions using global meteorological analyses such as from NCEP, but there are 
two reservations for doing this given the type of measurements associated with Cloud 
Slicing: (1) The analyses have a coarse resolution compared to the size of TOMS 
footprint measurements, and (2) the analyses must be accurate to within about 100 km to 
200 km in identifying the location of the windjets, otherwise a significant portion of bad 
data may be retained and good data dispensed with Cloud Slicing. 

Appendix B: Column Abundance and Air Mass in the Atmosphere. 



One can relate ozone and air mass in the atmosphere in a self-consistent way using the 
concept of column abundance and the Dobson Unit. The fundamental Dobson Unit 
definition applies to all constituents in the atmosphere, including the total atmosphere 
itself where both volume and mass mixing ratio are unitary (100 percent content). One 
Dobson Unit (DU) is equivalent to 10' 5 m vertical thickness of constituent under 
conditions of 

standard pressure (1013.25 hPa) and standard temperature (273. 16K) [e.g., Andrews et 
al, 1987, and references therein]. Under this definition, column amount (Q, in DU) 
between two pressure surfaces Pj and P 2 may be determined by integrating volume 
mixing ratio (X, in units ppmv) over pressure P (in units hPa) from 


Q = 0.79 f 2 XdP 


(Bl) 


The constant 0.79 in (Bl) was derived for upper tropospheric levels [Ziemke et al., 2001], 
but can be applied to any constituent in the troposphere or stratosphere with at most 
around 1% relative error. For the total atmosphere, Xis a constant (10 6 ppmv) in (Bl) 
which yields a simple relationship for computing the total air column lying between two 
arbitrary pressure levels P l and P 2 : Q = 7.9xl0 5 (P 2 -F ; ). For the global mean surface 
terrain pressure (about 990 hPa from NCEP re-analyses) this yields a global mean column 
amount Q of 7.8x10 s DU for the total atmosphere. 


Mass associated with column amount Q is determined by integrating local measurements 
of Q. over surface area element dS (units m 2 ): 

M = cr jodS 

5 (B2) 

It is well known that one Dobson Unit is equivalent to 2.69x1 0 20 molecules m' 2 [e.g., 
Ziemke et al., 2001, and references therein]. This number applies to any individual 
constituent in the atmosphere including the total atmosphere itself. In (B2), it follows 
that a is the surface mass column density (units kg m' 2 DU' 1 ) and is equivalent to 
2.69x10 2Q \xJN a , where p is mean molecular weight (48 for ozone, around 29 for the total 
atmosphere) and N A is Avogadro’s constant (6.022x1 0 26 particles kmol' 1 ). For the total 
atmosphere of the Earth, the column amount 7.8x1 0 8 DU from (Bl) applied to (B2) 
yields a total mass of around 5.2xl0 18 kg. In comparison, for ozone in the atmosphere the 
average total column is about 300 DU and the total mass of ozone is around 3.3xl0 12 kg. 

From (B2), a surface mass density associated with local measurements of Q is given by 
crQ. For NCEP atmospheric air mass this parameter becomes 7.9x \0 5 <j{P 2 -Pi). It 
follows that a proxy for lower stratospheric air mass can be obtained by replacing P 2 with 
tropopause pressure from NCEP re-analyses and choosing Pi as an arbitrary constant 
pressure lying above the tropopause. This investigation uses NCEP tropopause pressure 



as a first-order proxy of lower stratospheric air mass for cross-correlations of seasonal 
cycles with upper tropospheric ozone volume mixing ratio from Cloud Slicing (section 

We have incorporated NCEP re-analyses tropopause pressure which invokes the WMO 
2K km' 1 lapse-rate definition [WMO, 1957]. 
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Figure Captions. 

Figure 1 . Scatter plots of TOMS above-cloud column ozone (in Dobson Units) versus 
co-located THIR cloud-top pressure (in hPa) for a grid point located near Wallops, 
Virginia (38N, 75 W). The top three frames are three days in July 1984 and the bottom 
three frames are days in August 1984. The numbers shown in each scatter plot are the 
derived ozone volume mixing ratio and 2-sigma uncertainties (in parentheses) in units 
ppbv. 

Figure 2. Average number of Cloud Slicing measurements per month in this study (see 
Appendix). 

Figure 3. (a) and (b): Monthly ensemble time series for years 1979-1984 at several 
WOUDC ground station locations (indicated) comparing upper tropospheric ozone from 
Cloud Slicing with ozonesonde measurements. Units are ppbv. 

Figure 4. Seasonal cycles (from years 1979-1984) of monthly ensemble time series of 
upper tropospheric ozone volume mixing ratio (units ppbv) from several ground station 
locations (including stations plotted in Figure 3). Cloud Slicing measurements are solid, 

ozonesonde measurements are designated by circl e s . Cloud Sli cing measuremen t s 

include ±2 a statistical uncertainties. 

Figure 5. Seasonal averages (DJF, MAM, JJA, SON, indicated in each frame) of TOMS 
upper tropospheric ozone volume mixing ratio (units ppbv) between 60S and 60N derived 
from Cloud Slicing. Seasonal averages incorporated 1979-1984 Nimbus 7 TOMS/THIR 
co-located measurements of above-cloud column ozone and cloud-top pressure. 

Figure 6. Same as Figure 5 but for TOMS total column ozone (in Dobson Units). 

Figure 7. Latitude versus month zonally averaged seasonal cycles in upper tropospheric 
ozone volume mixing ratio (units ppbv) from Cloud Slicing. Seasonal cycles were 
derived from TOMS and THIR co-located measurements for 1979-1984. 

Figure 8. Similar to Figure 7, but for TOMS total column ozone in Dobson Units. 

Figure 9. Upper tropospheric zonal mean ozone volume mixing ratio for January 1979- 
December 1984. Units are ppbv. 



Figure 10. Comparisons of the annual (left) and semiannual (right) amplitudes in total 
column ozone (solid) and upper tropospheric ozone (dotted). Annual and semiannual 
amplitudes were derived by applying basic Fourier analysis to the 1979-1984 monthly 
ensemble time series data. Units for both annual and semiannual amplitudes are ppbv. 

Figure 1 1 . Latitude versus month zonally averaged seasonal cycles in lower stratospheric 
(46-100 hPa) mean ozone volume mixing ratio (units ppmv) from HALOE 
measurements. 

Figure 12. Latitude-dependent temporal correlations between upper tropospheric ozone 
and HALOE (dashed, stars), TOMS total column ozone (bold), and NCEP tropopause 
pressure (dotted). Correlation amplitudes exceeding 0.58 pass 95% confidence level (i.e., 
5% significance level). Seasonal cycles for upper tropospheric ozone, total ozone, and 
NCEP tropopause pressure were determined from 1979-1984 data. Seasonal cycles for 
HALOE lower stratospheric column ozone were evaluated for the 1993-2000 (i.e., post- 
Pinatubo) time period. 
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Popular Summary 

This study evaluates the spatial distributions and seasonal cycles in upper tropospheric 
ozone (pressure range 200-500 hPa) from low to high latitudes (60S to 60N) derived from 
the satellite retrieval method called “Cloud Slicing” . Cloud Slicing is a unique technique 
for determining ozone profile information in the troposphere by combining co-located mea- 
surements of cloud-top pressure and above-cloud column ozone. Measurements of Nimbus 
7 Total Ozone Mapping Spectrometer (TOMS) above-cloud column ozone and Nimbus 7 
Temperature Humidity Infrared Radiometer (THIR) cloud-top pressure during 1979-1984 
were incorporated. In the tropics, derived upper tropospheric ozone shows year-round 
enhancement in the Atlantic region and evidence of a possible semiannual variability. Up- 
per tropospheric ozone outside the tropics shows greatest abundance in winter and spring 
seasons in both hemispheres with largest seasonal and largest amounts in the NH. These 
characterises are similar to lower stratospheric ozone. Comparisons of upper tropospheric 
column ozone with both stratospheric ozone and a proxy of lower stratospheric air mass 
(i.e. , tropopause pressure) from National Centers for Environmental Prediction (NCEP) 
suggest that stratosphere-troposphere exchange (STE) may be a significant source for the 
seasonal variability of upper tropospheric ozone almost everywhere between 60S and 60N 
except in low latitudes around 10S to 25N where other sources (e.g., tropospheric transport, 
biomass burning, aerosol effects, lightning, etc.) may have a greater role. 


